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We introduce and solve a generalized model of 1-f-lD Lorentzian triangulations in which 
a certain subclass of outgrowths is allowed, the occurrence of these being governed by a 
coupling constant f3. Combining transfer matrix-, saddle point- and path integral tech- 
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is described by a ID quantum Calogero Hamiltonian. The coupling constant f3 survives 
the continuum limit and appears as a parameter of the Calogero potential. 
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1. Introduction 

Random triangulations are interesting both from a field theoretic and from a statistical 
mechanical point of view (for a review, see for instance 0-0). In particular, the so- 
called dynamical triangulations, also denoted in the following as Euclidean triangulations, 
provide us with a consistent lattice regularized version of 2D Euclidean quantum gravity, 
which in its scaling limit reproduces the continuum Liouville field theory. Moreover, the 
introduction of statistical mechanical degrees of freedom on dynamical triangulations also 
led to the discovery of new universality classes describing critical phenomena on fluctuating 
surfaces. 

Beside dynamical triangulations, another type of random triangulations has recently 
been introduced, known as 1+1D Lorentzian triangulations ||. What characterizes a 
Lorentzian triangulation is that it can be decomposed into slices by cutting along horizontal 
lines, as shown in Fig. 1. Each slice is made of an arbitrary sequence of triangles pointing 
up or down at random. Each triangle has two "time-like" edges within the slice and one 
"space-like" edge shared with a triangle in a neighboring slice. 
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Fig. 1: A Lorentzian triangulation (a) and its dual (b). 

Lorentzian triangulations were invented with the aim of constructing a model of quan- 
tum gravity where causality was built in at a fundamental level. The resulting gravity 
model was denoted as Lorentzian gravity @. This model has a well defined scaling limit 
which is different from that of Euclidean gravity. At present there does not exist any 
continuum formulation a la Liouville with such a causal structure. We note, however, that 
a seed of such a formulation might be found in reference [ffl . Lorentzian triangulations are 
also interesting in their own right as a new statistical model, allowing for the definition of 
a new class of lattice models. Lorentzian triangulations lie somewhere between regular lat- 
tices with defects and completely random lattices (i.e. dynamical triangulations) and it is 
possible that models based on this type of lattices will reveal yet other universality classes. 
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Fig. 2: A generalized Lorentzian triangulation in the dual picture. 

In this spirit we will study in the present paper a generalized model of Lorentzian triangu- 
lations where the time slice structure is preserved but where some outgrowths within the 
time slices are allowed. 

Such a generalized Lorentzian triangulation is displayed in Fig. 2 in the dual picture 
and is obtained by decorating an ordinary Lorentzian triangulation such as that of Fig.l (b) 
by adding arbitrary space-like "arches" connecting two points of the same time line. These 
arches do not intersect each other nor the vertical straight lines. For later convenience, the 
boundary condition imposed on the arch configurations is such that any number of arches 
can escape from the right hand side of a given slice and reappear on the left hand side, 
see Fig. 2. An elementary arch, in the language of triangulations, corresponds to a pair 
of connected triangles that do not propagate in time, as opposed to pairs of triangles dual 
to a vertical edge. In the terminology of quantum gravity such an object constitutes a 
particular type of (small) "baby-universe" . Keeping to the statistical mechanical language 
we will simply refer to such decorations as "outgrowths". Similarly, more involved arch 
configurations correspond to more involved outgrowths each living in a single time-slice. 
Clearly, the class of outgrowths we allow for constitutes only a small subset of the baby- 
universes present in the dynamically triangulated or Euclidean quantum gravity model. 
In order to get all baby-universes appearing in Euclidean gravity one would have to allow 
for loops connected in all possible ways to the already existing elements of the lattice via 
three- valent lattices. Such a model has been studied in ref. |J. It corresponds to the 
fully packed phase of the so-called 0(1) model on a random lattice. Note finally that the 
outgrowths we consider would not suffice to generate surfaces of arbitrary topology. The 
surfaces we consider have genus zero or one depending on boundary conditions. 

Our motivation for studying the particular type of lattices depicted in Fig. 2 is two- 
fold. First, in the context of 2D quantum gravity, it has been shown that one can view 



2 



Lorentzian quantum gravity as a renormalized version of Euclidean quantum gravity |J. 
Obviously, in going from Euclidean to Lorentzian triangulations baby-universes must be 
integrated over. A key point in the above renormalization procedure consists in a non- 
analytical redefinition of the "boundary cosmological constant" which is the parameter 
coupled to the length of the spatial slices, i.e. the length of the horizontal lines in Fig. 1 
or Fig. 2. In our model all outgrowths lie along the horizontal lines and one of our aims is 
to investigate whether the redefinition of the boundary cosmological constant can be due 
to these decorations alone. 

Secondly, from the statistical mechanical point of view, lattices of the type depicted 
in Fig. 2 are interesting because the system of arch configurations is known to be crit- 
ical |7| and by defining this system on Lorentzian triangulations one might be able to 
change the universality class of the geometrical system. So far one has not been able to 
study analytically the behavior of Lorentzian triangulations when critical matter fields are 
introduced. One has studied numerically systems consisting of one to eight Ising spins on 
Lorentzian triangulations [BJ . These studies show that the interaction between matter and 
geometry is much weaker than for the models based on dynamical triangulations. Analyti- 
cal investigations have been carried out for Lorentzian triangulations equipped with various 
dimer fields and for the case where higher curvature interactions are present 0. These 
models all turned out to belong to the same universality class as pure Lorentzian triangu- 
lations. Remarkably, the concept of integrability generalized to non-regular (Lorentzian) 
lattices played a crucial role in their exact solutions and led to the hope of extending the 
standard techniques of integrable lattice models to models on fluctuating lattices. 

The paper is organized as follows. In Sect. 2 we describe in more detail the generalized 
Lorentzian triangulations that we are going to study and write down their transfer matrix 
description. Next, in Sect. 3 we reformulate the model using a Schwinger-type integral 
representation leading to a continuous transfer matrix. The Schwinger language provides 
us with an alternative way of obtaining the exact solution in the case of pure Lorentzian 
triangulations. In Sect. 4 we consider the case where arches or outgrowths are present 
and show the existence of three different regimes of the model according to the value of a 
parameter f3 governing the density of outgrowths. In the interesting case of low density of 
outgrowths ((3 < 1), we define in Sect. 5 a sensible scaling limit in which our model gets 
identified with a well-known quantum mechanical system, namely that governed by one- 
dimensional Calogero Hamiltonian. Using this equivalence, we compute several interesting 
thermodynamic quantities. Finally, we discuss in Sect. 6 various aspects of our results. 
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2. The model 

Our generalized model is based on triangular lattices which in the dual picture look 
as in Fig. 2. In addition to a Boltzmann weight g per triangle, we introduce a weight h 
for each triangle which in the dual language is part of a vertical edge and a weight 9 for 
each triangle which in the dual language is part of an arch, see Fig. 3. We also introduce 
the parameter 

(2.1) 



governing the density of arch decorations. 
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Fig. 3: The weights associated with the different elements of a generalized 
Lorentzian triangulation. 

The partition function of our model over a time lapse t is obtained by summing over all 
such generalized triangulations with exactly t time slices and weighted as just explained. 
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Fig. 4: The transfer matrix G is decomposed into two parts: the matrix 
0W of the pure case, and the matrix G^ 2 -* implementing the upper and lower 
arch decorations. 
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As in the case of pure Lorentzian triangulations, we may describe our model using a 
transfer matrix approach. The transfer matrix describes the effect of adding one time-slice 
to our lattice. In the present case it is convenient to view the addition of one slice of 
the lattice as two consecutive steps. We first apply the transfer matrix of pure Lorentzian 
triangulations which corresponds to adding the portion of lattice below the dotted line 
and above the dashed line in Fig. 4, including a line of constant time. The transfer matrix 
element 0^ is indexed by two non-negative integers i and j which denote the number of 
half-edges going into the time line and that coming out of the line respectively and counts 
the (weighted) arrangements of these half-edges along the time line. We then complete the 
slice addition by applying a transfer matrix 0( 2 ) which corresponds to adding the portion 
of lattice above the dotted line and below the dashed one in Fig. 4 and accounts for the 

(2) 

pairing of half-edges into arches or into vertical edges. The transfer matrix element 0^ - 
is also indexed by the numbers of in and out-coming half edges. The total transfer matrix 
can therefore be written as a product of the two above transfer matrices 

e« = E e S )0 8 ) - ( 2 - 2 ) 

k 

In the case 9 = 0, h = 1 of pure Lorentzian triangulations, the matrix reduces to the 
identity. The new ingredient is therefore entirely contained in the definition of 0( 2 ). 

For convenience we will work mainly in the language of generating functions, by 
defining 

oo 

Q(x,y)= J2 (2-3) 

i,3=0 

and similarly for 0W and 0( 2 ). The composition law ( |2.2| ) then reads 

G(x,y)= I -^-e^ (x, -)0^ (u, y), (2.4) 

with the contour C encircling the origin. The generating functions can be easily derived, 
with the result 



e (1) (x,y) = - 

1 — qx 



gx - gy 



era,* v) = c ^ 2 ) ^V) 1 (L>5) 

K^yj i-e^x 2 c 2 {e 2 x 2 )i-e 2 y 2 c 2 {e 2 y 2 )i-h 2 xyC{e 2 x 2 )C{e 2 y 2 y 
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where C(z) is the generating function for Catalan numbers 



CM = ■ (2-6) 

The result for 0W is the direct consequence of the matrix element Q^) taking the value 
( t "i" , )fl'* + ' ? '> where we choose for convenience to include the weight g per triangle in Q^ 1 '. 
The derivation of G^ 2 -* is slightly more involved: we attach to each vertical edge the upper 
and lower arch configurations sitting immediately to its right (see Fig. 4). This results in 
an effective weight h 2 xyC(9 2 x 2 )C(9 2 y 2 ) per vertical edge. Indeed, C(z) is the generating 
function for arbitrary arch configurations with a weight z per arch. Summing over the 
number of vertical edges leads to the third factor in the second line of ( |2.5|) . The first 
two factors come from our particular choice of boundary conditions. Each is obtained 
by summing over the number of escaping arches each effectively weighted by a factor 
9 2 z 2 C 2 (9 2 z 2 ) , with z = x resp. y for upper, resp. lower escaping arches. We remark that 
as opposed to what is the case for 0^ and O^ 2 -*, the transfer matrix is not symmetric 
with respect to x and y.i 

Among other quantities of interest, we will compute the t-step partition function with 
periodic boundary conditions in the time direction 

Z(t) = Tr(e*). (2.7) 

We will also calculate correlation functions for the total number N(s) of triangles in the 
slice at time s, with the same t-step periodic boundary conditions. 

Apart from these we will be interested in calculating the t-step partition function with 
open boundaries in the t-direction Z(i,j, t) 



Z(i,j,t)=((0W0( 2 ))*0( 1 )).. ) (2. 



4 This asymmetry of Q(x,y) can be cured by, instead of defining G by ( |2 . 2| ) , setting 
= (0^ 1 - ) ) 1 / 2 0^ (O^) 1 / 2 . Since the transfer matrix of the pure model has been explicitly 
diagonalized Q it is straightforward to write down an expression for (9 (1) ) 1/2 . However, the full 
transfer matrix 0^ resulting from this definition has no longer a clear geometrical interpretation 
as being associated with i incoming edges and j outgoing ones. In appendix A we show how to 
calculate = (0^ 1 ^) 1/ ' 2 0^ 2 ^(0^ 1 ^) 1//2 . Here we shall take a line of action which does not require 
the knowledge of the exact transfer matrix. 
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or rather its associated generating function Z(x, y, t). With a slight abuse of language we 
will also refer to Z(x,y,t) or Z(i,j,t) as the loop-loop correlator. Notice that Z(i,j,t) is 
not simply the t-th power of the transfer matrix. It is, however, to our opinion the most 
natural definition of a £-step partition function with open boundaries in the t-direction. 
It describes triangulations for which the boundary conditions in the t-direction are as 
depicted in Fig. 2. There are i incoming vertical edges and j outgoing ones and no arches 
are allowed either between incoming or outgoing edges. 

For pure Lorentzian triangulations, two different strategies were used for calculating 
objects as given above. One strategy consists in first determining @ t (x, y) by writing down 
and solving a recursion relation in t for this quantity [[| . The other one consists in explicitly 
diagonalizing the transfer matrix 0^ 0. Neither method appears to be tractable for the 
present problem with general 9 and h. We will resort to yet another approach using an 
integral representation of the transfer matrix. 

3. Schwinger representation 

3.1. Transfer matrix 

In this section we will reformulate the discrete transfer matrix with indices i,j = 
0,1,2,... as a continuous transfer kernel by means of an integral representation, using 
continuous Schwinger parameters a, a' G [0, oo). This representation is inspired by the 
calculation in ref. || of the "watermelon" correlation function of the 0(n) model coupled 
to 2D Euclidean gravity. For definiteness, let us consider the partition function Z(t) (|2.7|). 
From fl2~lj) it follows that 

m = f[( /^^)n© (1) f— e^(, r ,«v), (3.i) 

AJ - \J C 2mw s 2lTTZ s J AJ ; \W r -i Z r J 

s=l x ' r=\ x ' 

where the contour is a small circle around the origin and where we have imposed cyclic 
boundary conditions on the w's, namely wq = wt- 

Whereas in the pure case the integrand in ( |3.1| ) would have singularities only in the 
form of poles, in the present case the integrand in addition has singularities in the form of 
cuts. We can, however, transform the cuts into poles by the following change of variables 



lj r = 9w r C(6 2 w 2 r ), Pr = 6z r C(6 2 z 2 ). 
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(3.2) 



This change of variables maps the circle C into a closed ellipse-like curve 8 still encircling 
the origin. Using the quadratic equation satisfied by the Catalan generating function 



x(C(x)) 2 = C(x) - 1, (3.3) 



the relation (|3.2|) is easily inverted into 



w r = — , z r = — 3— . (3.4) 



From there it follows that 

dw r 6(1 — /J%)dfi r 



w 2 fi 2 . 



(3.5) 



and a similar relation holds between dz r and dp r . Inserting ( [3.2|) and ( |3.5| ) in (|3.1| ) we get 

m = n (i 2^2^:) n i- 1 ^i-,% r+ ^ 1 + ^_ 1+ ^)' (3 - 6) 



where (3 = 6 /h as in ( |2.1| ). We may next introduce a Schwinger representation for the 
second factor in each of the products in the integrand. This gives 

The integrals over the p r may now be performed, by noticing that each of them simply 
picks the residue of the integrand at p r = This gives 

m - n (i ^if — ) (^) V— ±> 



(3.8) 



where the I^s are the modified Bessel functions of the first kind, i.e. Ifc(2a;) = 
S P >o x 2p+k /(p\(p + k)\) and where we have imposed periodic boundary conditions on 
the ct's, i.e. at+i = ot\. 

Similarly, to carry out the integral over \i r we simply have to pick the residue of the 
integrand at \i r = and the expression finally reduces to 

Z ® = II / da * e ~ aa IT M90a r ,gea r+1 ), (3.9) 
s=1 \Jo J r=1 



where we have defined the transition function 



M^V) = ^4(2x)/ fc (2y)//3 2fc . (3.10) 



fc>0 



Note that as expected the partition function depends on g, h and 9 only via g9 and (3 = 9/h. 
Eqn. (|3.9|) leads to the definition of a continuous symmetric transfer kernel 



G/3 !g o(a,a') = e <f)p(g9a,g9a'), (3.11) 

to be integrated with the flat measure on the parameters a > 0. To make contact with the 
original discrete formulation of the problem, let us show how to obtain statistical properties 
involving the observable N(s) = i + j counting the total number of triangles in the slice 
at time s (i and j denote respectively the number of in- and out-coming half edges at the 
time-line s), from the statistical properties of the variable a s . Let us define a more general 
local observable E s (z) through 

a n (s) = N(s)(N(s) - l)...(N(s) -n + 1), 

~ z n (3.12) 

n=0 

The correlations of E at various times s\, s^. read simply 

k k Zm 

( E Sm (^ m ))e = ( I] "-T^ 1 -)^ (3.13) 

m=l m=l 

where the subscripts G, G refer to the framework (discrete, continuum) in which the 
correlation is evaluated. This is readily proved by noticing that the desired correlator 
corresponds to substitutions g — > g(l + z m ) in the formula for the partition function 
( p.9|) within the slice s m . The latter can then be absorbed into a change of variables 
a Sm — > a Sm /(l + Zm). This will be extensively used in Sect. 5 below. 

3.2. Example 1: pure Lorentzian triangulations revisited 

As a preliminary exercise, it is instructive to re-derive the exact results of pure 
Lorentzian triangulations in this new continuum language. Let us diagonalize the con- 
tinuum transfer matrix of pure Lorentzian triangulations (with 9 = and h = 1) 

G (0) (a,a') = lim G^ g p{a,a') = I {2gVao'), (3.14) 
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where we have used that Ik(2gj3x)/ (3 k — > (gx) k /kl in ( 3. ICQ . It turns out that the eigen 



functions of are simply related to the Laguerre polynomials L n (x), orthogonal wrt. 
the measure e~ x dx over the positive reals. These polynomials read 

' \m I ml 

m=0 V 7 

The diagonalization follows from the following quadratic formula 

J2e--L n (x)e-U n (y)u^ = ^e'^^I, ( ) , (3.16) 

z — ' 1 — u z \ 1 — u z I 

n=0 V 7 

valid for u < 1. Hence upon setting g = l/(q + 1/q), and defining a = x(l + u 2 )/{l — u 2 ), 
a' = y(l + u 2 )/{l — u 2 ), with u = q, eqn. (|3.16|) translates into 

00 

G^(a,a') = ^Vi 0) («)^ 0) («')An, (3.17) 

n=0 

in which ifjn\a) is the normalized eigenfunction for the eigenvalue A n , respectively reading 



g 2n+l (3.18) 



A, 



5 

2^ 



9 = gC(g 
with C as in (gg). 

This leads immediately to the pure Lorentzian triangulation partition function on a 
time cylinder with t steps 



n=0 



Similarly, we get the one-point correlation of the observable T, s (z) as in ( |3.13| ) by using 
the formula (|3.16| ) with x = y = a(l — q 2 )/(l + q 2 ) and u = q tm . 

(S.(z)> e = ^(e— ) G( o) = ^ (3.20) 



where 



I(z) = — dae a y—* ^^>I 2a 9 2t \ ■ 3.21 
1 + 2 Jo V 1 -<? ! + <? / 
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Performing explicitly the integral, we find 



1 — g 1 1 — g- 4 y 1+3 1—9 

The expressions ( |3.19| ) and (|3.22| ) have a singularity at q = 1, i.e. g = 1/2. In order 
to obtain a continuum theory, we therefore write 

2g = l- ^a 2 A, (3.23) 

where a 2 is a scaling parameter with the dimension of area and where A is the renormalized 
fugacity per triangle. Using this scaling for g we have 

q = g C(g 2 )~ e - aVX , (3.24) 

and we see that in order to obtain a finite expression for Z^ (t) in the continuum limit we 
must set 

£=-. (3.25) 
a 

This gives the following continuum partition function 

4°) = lim ±z(°\t) = — ^. (3.26) 

T a^O 2* W 1 _ e -2v / AT V ' 

Similarly, we must take 

z = aVIZ and N =—, (3.27) 

a 

where the factor a/A is simply for convenience (Z is dimensionless) , and substitute this 
into ( |3.22| ) to get the scaled average 

(e^ 2 ) = 1 (3.28) 



1 - 2Zcotanh( v / AT) + Z 2 

generating the moments (J\f k ) of the (rescaled) number of triangles per time slice. These 
moments read 

(X k ) = ^ r P fc (cotanh(v / AT)) = -^ Lm=oU _ , (3.29) 

where Pk(x) denotes the k-th Legendre polynomial. 
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3.3. Example 2: the case (3 = oo 

For (3 — > oo the continuum transfer kernel factorizes as follows 

Gocgeia, a') = e -'/ (2^a) x I (2g6a'). (3.30) 

The partition function therefore factorizes too and reads 

= (Z^l))*, 

2 „ v2 . (3-31) 



poo 2 

Z^(l) = / dae- a I (2g6a) 2 = -K{{Ag9) 
Jo n 



where K{x) is the complete elliptic integral of the first kind. Note that Z(°°\l) is simply 
the partition function of two (upper and lower) interlocking arch systems connecting points 
by pairs along a single time-line. The function K(x) is singular at x = 1. In order to define 
a continuum theory we would therefore in this case set 

Ag9 = 1 - a 2 A. (3.32) 

Then we get 

Z(°°\l) = -K((l - a 2 A) 2 ) = --Log(a 2 A) + (9(1), (3.33) 

7T 7T 

and the leading singular behavior of (t) is therefore 

Z^\t) ~ ^-iLog(a 2 A)^) . (3.34) 
Thus, in this case it is not possible to define a sensible continuum time variable. 

4. Scaling regimes 

It is obvious that the possible singularities of Z(t) must come from large Schwinger 
parameters, a r , since the Bessel functions /fc(2x) are polynomial for small x and behave 
for large x as 

We will therefore be interested in determining the asymptotic behavior of the transfer 
kernel Gp^ g e{a, a') as a, a' — > oo. It turns out that one has three different scaling regimes 
corresponding respectively to (3 < 1, f3 = 1 and /? > 1. 
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4-1. Three scaling regimes 

When (3 > 1, the transition function (f>/3(x,y) in ( |3.11] ) is expressed as an absolutely 
convergent series (]3.1U|) , that behaves for large x and y as 



e 2(x+y) 1 ^2 e 2(x+y) 

In order to determine the asymptotic behavior of <f>/3(x, y) for the remaining values of (3 we 
make use of the generating function for modified Bessel functions of the first kind 

e*< t+ *> = I (2x) + J2(t h + ^)h(2x). (4.3) 

A;>1 

From this relation it follows that 

r dt eX{t+ 1 )+ y(pH+^- t ) = ( y) + ( } _ /o(2x)/o(2y ) 

= i (2^{x~TWyW+ym) ■ 

For /3 = 1 we immediately get from ( |4.4| ) 

01 (x, y) = i (I (2(x + y)) + I (2x)I (2y)), (4.5) 
which behaves for large x and y as 

fcfay) ~ ——==. (4.6) 
4 v /tt(x + y) 

Finally, when (3 < 1, we write ( [4.4|) as 

<t>f,(x,y)=Io (2y/(x + ^y)(x + y/^ + I (2x) I (2y) - 0i (x, y), (4.7) 

and notice that since <fii{x, y) is now the absolutely convergent series, the last two terms 
in ( P~7 ) behave like e 2 ( x+y - > . Using moreover that 

y/(x + [3 2 y)(x + y/[32)>x + y, (4.8) 

for all real (3^1, the large x, y asymptotics of 4>p(x,y) is entirely governed by the first 
term on the rhs. of ( |4.7|) . This means that for (3 < 1 we have the following asymptotic 
behavior of the transition function as x, y — > oo 

(^(x,y) I 1 + +... . (4.9) 

2^/it^(x + [3 2 y)(x + y/f3 2 ) \ 16^ (x + /5 2 y)(x + y//? 2 ) J 
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To summarize, we have found the following dominant behavior for the transfer kernel 

G/3, g 9(a,a') ~ e - s ^^^') as a? a ' _> qq. 



P>1 : Sp, a9 {a,a') = ~(a + a')(l - 4g9), 



1 



/3<1 : S/3,go(at,a') = -(a + a') -2gd x j{a + — )(a + P 2 a'), 



(4.10) 



where the cases (3 > 1 and /3 = 1 differ by sub-leading corrections. A crucial difference 
between (3 < 1 and f3 > 1 is that in the latter case the partition function is factorized 
(at leading order), while in the former case it will be dominated by correlated en's all of 
the same order. This in turn reflects the natural property that for low arch densities, 
successive time slices have lengths of the same order, while the presence of many arches 
allows for decorrelated lengths. 

The following three sections are devoted to the three cases (3 = 1, f3 > 1 and f3 < 1 
for finite t. In all cases, we will need the large-a asymptotics of the product of t transfer 
kernels which may be expressed as e~ Sl3 x Up, with an action Sp({a s }), reading 

t 

Sp({ a s}) = ^2 S P:go{ a s, a s +i), (4-11) 

s=l 

where ctt+i = a ii an d a non-exponential factor Up({ot 3 }), which will be detailed below. 

4.2. The case (3 = 1 

The action S\({a s }) and the pre-factor Ui({a s }) read respectively 

t 

S 1 ({a s }) = (l-4 9 e)J2^s, (4.12) 

i=l 

and 

Ui({a s }) = \ fl , = ( 4 - 13 ) 



From the action ( |4.12 ), we immediately find the singularity of Z(t) to be at 

Ag9 = 1. (4.14) 

Furthermore, writing 4g0 = 1 — a 2 A and performing the change of variables a s = f3 s /a 2 in 
the integral ( |3.9| ) with j3 = 1, we get 

a*J Q ^2^ + ^+1 
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when a — > 0. Note that the multiple integral converges and gives a non-trivial dependence 
on t. 

This calculation should be compared to that of computing the singularity structure 
of the t-slice "watermelon" correlator of the 0(n) model coupled to ordinary 2D Euclidean 
gravity. Indeed, viewing the slices of watermelon as time-slice separators, the correlator 
becomes in the limit n — > equivalent to a periodic time Lorentzian partition function 
with h = and with slightly different boundary conditions (no escaping arches). In both 
cases the singularity of the partition function is factorized, leading in the case of ref. 



to a different behavior Z ~ a , in agreement with the KPZ scaling [10 



4.3. The case (3 > 1 



In this case, the action is the same as in the /3 = 1 case Sp = Si (|4.12 ), while the 
pre- factor U/3({a s }) reads 

u >«*» = iw^w)' & (416) 

The critical point is therefore the same as in the case (3 = 1, namely 

4</0 = l, (4.17) 
and setting again 4g9 = 1 — a 2 A, the singularity of Z(t) is factorized with the result 



1 (3 



2 



t 



which, up to a proportionality factor, is identical to the behavior at (3 = oo (|3.34|) . 

We may infer that this behavior is characteristic of the (3 > 1 regime, and breaks 
down at /? = 1 as ( [118| ) indicates. Once again, for (3 > 1, it is not possible to introduce a 
sensible continuum time variable. 

4-4- The case (3 < 1 

In this last case, the dominant action reads 

t t 



SpUpt.}) = - 2g0j2\l( a * + ^s+i)(a s + (4.19) 



s=l s=l 
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and the sub-leading factor is 



1 + / 

16g9j(a 8 + (3*a s+1 )(a s + ^) 



(4.20) 



where in both expressions we assume periodic boundary conditions for the a's, i.e. ott+i = 
ol\. Let us now turn to studying the singular behavior of Z(t). To obtain the critical 
point, we notice that the dominant action ( |4.19| ) is homogeneous of degree one in the 
a's, thus in order for the partition function to be defined, all configurations {a} must 
satisfy Sp({a}) > 0. The critical point corresponds to the existence of configurations {a*} 
minimizing the action, while satisfying Sp({a*}) = 0. These can be obtained from the 
saddle-point equations 

!,f) f(3 2 x 8 + — + x s _i + ) = 0, s = l,2,...,£, (4.21) 



ft \ ~ 

p y X s Xg — i 

where 



x s = ( id 1 — ~ ) ' Xs +t = Xs and CXs + t = as - (4.22) 



fl 2 a s + a s+1 

The only periodic solution a s to the above corresponds to x s = const. = 1 for all s, 
implying that all the a's must be equal, and moreover this forces g9 to take the critical 
value 

2 9 o(p+^\ =1. (4.23) 

As expected the action at the saddle-point vanishes automatically. 

To capture the singularity of Z(t), let us explore the vicinity of the critical point 
( (4.231 ). In analogy with the pure triangulation case, we set 

2g9 ((3+1^=1- ia 2 A, (4.24) 

where a 2 is a small parameter with the dimension of area. The homogeneous nature of the 
saddle-point above, where all the a's are equal, suggests that we keep, say a\, unchanged, 
and that we perform the change of variables a s = a\ + £ s ^cti for s — 2, 3, t. Then the 
sub-leading factor (|4.20| ) together with the integration measure behave for large ol\ as 

fl da.Up{{<*.}) - ~P= II (4-25) 

s=l S=Z 
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while the action reads 



s ( a i, {£}) = ta i + 



8 

S=l 

t 



where we have set £1 = = 0. To evaluate Z(t), we must first integrate over the £'s 
(along the real line) and then over oti > 0. This results in 

l-(3 2 fl+P 2xt ~' 



2aVAt(l + f3 2 ) V 1 - P 

where the £ integral has produced the term t -1 / 2 = det(A) -1 / 2 , where A is the (t — 1) x 
(t — 1) matrix of the discrete Laplacian, namely A rjS = 25 r>s — 8 r}S +i — 5 r ,s— ij truncated 
to 1 < r, s < £ — 1. Note that in (|4.27|) we have displayed the trivial entropic factor 
(2(1 +/3 2 )/(l -(3 2 )Y, namely a weight 2 for each application of (1 > and (1 + (3 2 ) / (1 - (3 2 ) 
for each application of G*- 2 -*. These factors could of course be absorbed in a redefinition 

e(D - eW/2 and e( 2 ) - e( 2 )/ (i±f^). 



5. The scaling limit and the Calogero Hamiltonian 



In order to obtain a finite result when a — > for the partition function (|4.27|) (after 



stripping it from the trivial entropic factor) , we must let simultaneously t scale as l/o. We 
therefore introduce a continuum time variable T through 

T 

£=-. (5.1) 
a 

We note that ( |5.1| ) is the same scaling of the time variable as in the pure case (cf. 
eqn. (|3.25|) ). This means in particular that the fractal dimension, dn of our surfaces 
for (3 < 1 is also the same as in the pure Lorentzian case, i.e. dn = 2. 
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5.1. Partition function 

Let us now compute the partition function of the (3 < 1 model in the vicinity of the 
critical point using Q4.24| ) and (|5.1| ) and letting a — > 0. We know already that the integral 
( pj.9| ) should be dominated in this limit by large ct's that are close to one another. This 
suggests that we directly perform the following change of variables in the integral 

2 

a s = ^, ip s > 0. (5.2) 
a 

Inserting the change of variables in the action ( |4.19|) leads to 



G , \ !r 2 (1 - |a 2 A) A 2 , ip 2 s+1 -ip 2 (p 2 s+1 - tp 2 

s=l x ' s=l r * 



2 i 



(5.3) 



2 

s=l 

where both terms are of order a, as <£>s_|_i — v^s is of order a. The contribution from the 
pre-factor Up({a s }) ( |4.2Q| ) and the measure reads 

n ^({a.}) = n i/J^ 1 + A) - (n -J^.) et EL < b - 4 > 

s=l s=l ^ s \s=l ^ / 

We notice that this latter relation has the effect of adding an extra term to the action 
( |5.3|) . Next, we introduce a discrete function by ip(u = s/t) = (p s and we assume that 
for large t (small a) this function becomes a smooth function of the continuous variable 
u G [0,1]. This allows us to Taylor-expand <^ s +i around (p s as (p(u + a/T) = <p(u) + 
^<p'(u) + 0(a 2 ). Furthermore, we can replace sums by integrals, J2 S ~^ 7T Jo 1 ^ Um Finally, 
it proves convenient to define the continuum partition function Zt by 

Z T = lim f-^lTl Z(t). (5.5) 



>o \2 1 + /3 2 

Inserting everything and rescaling <p — > ( ) <P, u ~ ^ we get 

£Ve 2 j o ^ w w I^T) 3 ^, (5.6) 



fT 

V(0)=y(T) 

where we have identified the functional measure as 

t 



Vy = lim T]-^=, ( 5 - 7 ) 

a s=l 
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and where the integration extends over real positive fields <p(u) defined for u G [0, T] and 
obeying (p(0) = <p(T). The parameters to and A are given by 

Using the Feynman-Kac formula we see that computing Zt amounts to solving a one- 
dimensional (imaginary time) quantum mechanical system on the real half-line, namely 

Z T = Tr(e" T ^), (5.9) 

with a Hamiltonian 

"=-\i?+r 2 ^-w' * >0 - <5 - 10) 

This is nothing but the celebrated Calogero Hamiltonian for one particle. We thus arrive at 
the result that 1 + 1-dimensional Lorentzian triangulations (in the low arch density regime 
(3 < 1) are equivalent in the continuum limit to a one-dimensional Calogero model. This 
holds in particular for the pure case. A Hamiltonian describing the continuum limit of 
pure Lorentzian triangulations, i.e. /3 — ^ was derived in ref. using a strategy different 
from the one employed here. The Hamiltonian of ref. |§ is expressed in terms of a variable 
L which has the interpretation of the continuum counterpart of the length i of a line of 
constant time. As we shall see later the field if 2 has the interpretation of the continuum 
counterpart of the total number of triangles in a given time slice % + j, also equal to the 
total length of the two adjacent constant time lines, translating into 2L in the continuum 
limit. Due to the use of different boundary conditions in ref. || and here one cannot 
immediately compare the two Hamiltonians. In appendix A we show how to recover the 
Calogero Hamiltonian for f3 — > using the approach of ref. [Q. As expected a change 
of variables 2L — > ip 2 is involved. We also explain why the approach of ref. does not 
extend to f3 ^ 0. 

Note that the parameter A of our Hamiltonian ( |5.10|) satisfies < A < 1 and that A = 
1 iff P = 0. It is well-known [|TT] that A = 1 is a limiting case for the quantum mechanical 
system ( |5.1U[ ). For A > 1 the operator H is no longer self-adjoint. Our generalized 
Lorentzian triangulation model thus always leads to a physically acceptable Hamiltonian. 
For [3 < 1 the Hamiltonian H is readily diagonalized as follows [0|. Introducing the 
confluent hypergeometric function 

. \ a(a + l)...(a + k — l) k ,~ 
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the eigenvectors of H take the form 

M<P) =A n e-^ 2 <p»-?F(-n,fi,u;<p 2 ), (5.12) 
where A„ is a normalization factor 



^-^nnlm»r (5 - 13) 

ensuring the orthonormality of the eigenbasis ijj n (ip), and the parameter \i may take two 
values, related to the potential strength A through fx = 1± hyl — A. However, the smaller 
of these two values must be discarded as unphysical leading to 

1 + /3 2 



^ = 1 + T^- ( 5 - 14 ) 



Note that the confluent hypergeometric series F is truncated by the value a = —n to a 
polynomial of degree n, generalizing the Laguerre polynomial L n , namely 

The corresponding eigenvalues read 

E n =uj(2n + /j) =u{2n + l + — ^—,), n = 0,1,2,... (5.16) 

1 + p z 

This leads to the scaled partition function 



e 



n>0 



For /3 — > the result ( |5.17| ) reduces to the pure Lorentzian triangulation result ( |3.26|) . 
On the level of eigenvalues and eigenfunctions we have lim a ^o(fi r ^n) t = lim^-a e~ EnT 
and lim a ^o ^n\ot)yfda = lim M ^i ijj n (tp)y/dip, respectively obtained from eqns. (|3.18|) and 
( |5.12| ) with to — > y/X, and after setting q = e~ a ^, t = — and a = as before. 

Furthermore, we see that as f3 — > 1 the exponents in ( |5.17| ) become singular, signaling 
another type of scaling coming into play (c.f. Sect. 4.2 and 4.3 above). It is worth stressing 
that the parameter (3 cannot be absorbed by a redefinition of the renormalized triangle 
fugacity A or the continuum time variable, and is itself a genuine continuum parameter. 
This will prove even more visible when we compute correlation functions. Such a situation 
where a coupling constant introduced at the discrete level survives in an unrenormalized 
form in the continuum limit is quite unusual. 
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5.2. The loop-loop correlator 

The Laplace transform of the loop-loop correlator Z(x, y, t) can be written as (cf. 
eqn. (£§)) 

Z(x,y,t)= Z(i,j,t)x l y J 

i,j>0 

t 

=n 



8=1 



dw s dz s 
q 2mw s 2mz t 
t-i 



x 



V W lJ \r=l ^ r W r+1 J ) \z t J 

(5.18) 

Repeating the computation of Sect. 3, namely performing the change of variables «v, z r — ► 
UriPr (cf. eqns. ( |3~^ ) and (|3.4|) ), representing the O^-factors as Schwinger-type integrals, 
and performing the explicit contour integrals over the p r and fj, r , it appears that a: and y 
are simply spectators throughout the computation, and one finds 

t / poo \ t— 1 

Z(x,y,t) = 11 ( / rfase " Qs ) e 9(aoa:+at?/) IJ^^«r,^a r+ i), (5.19) 



s=0 v,yu 7 r=0 

where the transition function <f>p(x,y) has been defined in ( 3 . 1 0|) . Apart from the explicit 
exponential dependence on x and y, the loop-loop correlator differs from the partition 
function ( |3.9| ) by involving one more Schwinger parameter. Note also that in this case we 
do not impose periodic boundary conditions on the en's. 

We will now study the behavior of the loop-loop correlator in the scaling limit defined 
by (|4.24|) and (|5.1| ). Here we must also assume that x and y are close to their critical 
values x c = y c = (corresponding to the poles of the O^-factors in (|5.18| )). It turns out 
that the correct scaling ansatz reads 

x =L(l-aX), y=-L(l- a Y). (5.20) 
To calculate Z(x,y,t) in the scaling limit we use the same strategy as in Sect. 5.1. We 

2 

perform the change of variables a s = s = 0, 1, . . . , £, introduce a discrete function by 
ip(u = s/t) = f s and assume that as t — ► oo this function becomes a smooth function 
of a continuous variable u G [0, 1]. However, in the present case we must deal with fixed 
as opposed to periodic boundary conditions for <p and special care has to be taken in the 
treatment of the boundary terms. Writing the integrand in ( |5.19| ) for large a s as e~ s<3 Up 
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in the same way as in Sect. 3, we get, when inserting our change of variables, the following 
contribution from the measure and the pre-factor 

(5.21) 

Inserting the various scaling relations in Sp({a s }) on sees that the scaling ansatz (|5.2Q|) is 
exactly what is needed to make divergent boundary terms cancel. Furthermore it becomes 
natural to define a continuum loop-loop correlator Zt(X,Y) by the following recipe 

Z T (X,Y) = \im a -J^(l^) Z(x,y,t). (5.22) 



Here we choose to scale away the same entropic factor as in the case of partition function, 
namely a factor 2 for each occurrence of the transfer matrix O 1 - 1 ) and a factor ( fz^p j 
for each occurrence of the transfer matrix O^ 2 ). Collecting all the terms and rescaling 
V 9 — *■ ( l+ps ) '■P an d u ~ *■ Tu as before, we find that the continuum loop- loop correlator can 
be expressed as the quantum mechanical propagator 

Z T (X,Y) = {^rpj J o °° ^dvo^d^e-MX^T^fale-TA^ ( 5 . 23 ) 

where the Hamiltonian H is defined as in (|5.10|) , with to and A as in (|5.8| ), and where we 
have set 

To evaluate Zt(X,Y) we simply have to insert two decompositions of the identity as a 
sum over projectors on the eigenspaces of H, i.e. I = ^2 n>0 \n)(n\ and use (n\(p) = ip n ((p), 
where i^ n { i p) is the n-th normalized eigenfunction of H, defined in ( |5.12| ). This yields 

Zt(X,Y) = (\^) 2 f:e- TE -G n (X)G n (n 

V P 7 «=o (5.25) 



G n {Z)= I ^d<pe-^ 2z ij n (<p) 



with E n as in ( |5.16| ). The function G n (Z) is readily determined using the expression ( |5.12j ) 

for ip n (ip) 

poo 

G n (Z)=A n / ^ e -^ 2 (-+^)F(-n,/i,^ 2 ), (5.26) 
Jo 
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where \x = 1 + \+[p as before. Inserting the polynomial expression for the truncated 
confluent hypergeometric function ( |5.15|) , eqn. ( |5.26|) can be integrated term by term with 



the result 



G n {Z) — l —-—^ J - £ 



2 \u + ZJ £f \kj //(// + l)...(n + k - 1) V u + Z 
An + 1 f 2 2^ , 



(5.27) 



2 v 2 y V^ + z / 2 w + 

Here we have recognized the ordinary hypergeometric function 

w , \ , _ . ^ fc a(a + l)...(a + fc-lW&+l)...(& + fc-l) 

F(a, 6, c, ,) S2 F l( a, 6, c, ,) = g>' -J i. (5 .28) 

Thus, the continuum loop-loop correlator reads 



z T (x,y) 



1 + /3 2 V 1 /r(^±i)\ / 4u; 



2 _ 



,2 \ — 



i-(3 2 J 2oo y r(/x) y V(^ + x)(^ + y). 

4- r(n + l) 1 ' 2 ,M 'u; + X j 1 ' 2 ' M ' 



(5.29) 



This expression may further be simplified by use of the following quadratic relation due to 
Meixner 



»=° (5.30) 

-F{a,b,fi, 



[l + s(l-z)) a (l + s(l-w)) b v ' (l + a(l-z))(l + a(l-w)) 



Indeed, applying (|5.30|) with s = — e 2Tw , a = b = , z = 2u/(u + X), and w = 
2u/(u + y), the continuum loop- loop correlator finally reads 

(^'^e^^^ ( , 31) 



where 

a; a; 



(X sinh(wT) + u cosh(wT)) (Y sinh(u;T) + u cosh(u/T)) 
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(5.32) 



Let us consider the pure Lorentzian triangulation limit (5 — > 0, i.e. (x — > 1 and u; — > 
\/A. This limit involves the hypergeometric function .F(l,l,l,i>) = 1/(1 — u), and the 
corresponding loop-loop propagator reads 



4°)(x,y) = sinh(wT) -± 



CJ 1 — 1) 



^ V " X _ (5.33) 
A v ; 



(XF + A) sinh( v / AT) + >/A(JC + Y) cosh( v / AT) ' 

This coincides precisely with eqn. (2.29) of ref. 0. 

Let us now translate the result ( |5.31|) into the language of conjugate rescaled loop 
lengths by setting 

Lq = ai, Lt = aj, T = at, (5.34) 
in terms of which the rescaled loop-loop correlator reads 

Z T (L ,L T ) = lim \^ (iZ^)' Z(i,j,t). (5.35) 

This is obtained by taking the inverse Laplace transform of ( |5.31|) over X and Y . In 
practice, it is more convenient to read it directly off the original expression ( |5.23j ) upon 
performing the change of variables 

Lo = ) , = 2 ^ \T=p ) ' (5 - 36) 

relating (p\ resp. to the numbers of triangles pointing up resp. down in the time slice 
u = resp. u = T. This yields 



Z T (L , L T ) 





2 (<Pt 


e -TH 






y/<PQ<PT 



(5.37) 



2 



To evaluate the above heat kernel, we can use the following quadratic relation satisfied by 
the generalized Laguerre polynomials introduced in Sect. 5.1, and generalizing the pure 
case ( |3.16| ), namely |D| 



E ^iZrWn-*, ft *)H-n, ft y)z" = ^-e-^i,yz)^ (5.38) 

Z — ' Till (/Uj^ 1 — 2 1 — Z 



n>0 
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where I„-i(2x) = J2 p>0 x^ 1 ^ /(p\T(p + //)) is the modified Bessel function. This equa- 
tion is immediately rephrased in terms of the eigenfunctions ijj n of H (|5.12|) as 



z^HMMvt) = 2u^^^e-^^^I,_ 1 U^^) . (5.39) 

n>0 Z \ Z J 

Taking z = e~ 2wT we deduce the expression for the heat kernel of the Calogero Hamiltonian 
H, leading finally to 

1 ~ P 2 \ 2 ^ -(L +LT)(i^-) 2 ^cota n h(a;T), 



Zt(L °' Lt) ~ ) imh(^T) 6 

f fl-P 2 V - X (5 ' 40) 

with cj as in ( |5.8| ) and fi — 1 = (3/(1 + /3 2 ). This result coincides exactly with eqn. (2.32) 

of i- 

5. 5. Correlation functions on a time cylinder 



In this section, we will compute for f3 < I the general correlations (|3.13| ) of the 
observables H Sm (z m ) (|3. 12 ) , in the continuum limit where a — > 0, with 



Sm = — , z m = auZ m and a Sm = . (5-41) 

a a 



With these substitutions, eqn. ( |3.13|) becomes 



k k 

2 



m=l m=l 



where the rescaled number of triangles in the time-slice u G [0, T] reads J\f(u) = 
a { i+p i Y 'N(u/a). Apart from its instrumentality, eqn. ( |5.42|) yields the interpretation 
of the field (p(u): the square of f{u) is identified with the rescaled number Af(u). This is 
very reminiscent of the field theoretical representation of polymers, where ip 2 is identified 
with the polymer density. For the pure case, this should not come as a surprise as we have 
shown in that Lorentzian triangulations can be mapped onto random walks. 
Let us first compute the one-point average of e u ^ 2 : 



f,{Z)^{e^ z ) = ^- d^^) 2 e- T ^e z ^\ (5.43) 

ZT ^ J ° n>0 
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with the cylinder partition function Z?^ = e~ wTtI / (1 — e~ 2ujT ), indexed by fi = l + + 
(3 2 ) for convenience. We use again the above quadratic equation (|5.39|) with z = e~ 2wT 
and v?o = ¥t = <P, and we perform the change of variables u = unp 2 to get 

U(Z) = e-^-i) r due -(cot^T)-z )ui f u \ 
Jo \smh(uT)J 

This last integral is a particular case of the generic Laplace transform of the modified 
Bessel |13|| : 



r 



/ due—I^su) = i 1 \ (1 - yfl - {8/r)*y-\ (5.45) 
Jo \/l - Is/rV 



and we finally get 



/ cotanh(ujT)-Z-^/l-2cotanh(u;T)Z+Z 2 \ ' 
I cotanh(wT)-! J 



/mO 2 ) = ~ > ( 5 - 46 ) 

M v/l-2cotanh(wT)Z + Z 2 



generalizing the pure case result ( |3.28| ) corresponding to \x = 1 and cj = vA. This yields 
alternatively the moments (Af k ) of the rescaled number of triangles per time-slice 

,| ,| y^fc /fc+^-lN rk+l—fj,\ —2mu)T 

= ^^"-"-"'(cotanh^T)) = ^'"=° , (5 . 47) 



where Pjf '(x) denotes the /c-th Jacobi polynomial |i~3| . 



The formula (|5.46| ) generalizes nicely to the case of the /c-point function ( |5.42| ), and 
we leave the details of its derivation to appendix B below. Defining the time intervals 

t m = u m+ i - u m , (5.48) 

for m = 1, 2, k — 1 and tk = T — Xa<m<fc-i ^ ne resu lt reads 



m=l 



where 



gii {s) = S mh(uT)e» T ^- l h 2 - ^ 77==^ = A* (cotanh(o;T) - - ) , (5.50) 
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with as in ( |5.46| ), and where P is the following polynomial of Z±, Z^ 



P(Zi, ...,Zk\h,...,t k ) = 2cosh(u;T) 

k r (5 51) 

+ E(~ 1 ) r E 2 mi Z m2 ...Z mfc n2sinh(u;( *-))> 

r=l Kmi <rn,2 < . . . <m r <fc J = l Tij <s<fHj^i — 1 

with suitable boundary conditions on the indices, namely: m r+ i = mi + k and t s = t s _k 
for s > k. Note the remarkable fact that all the /c-point correlators are expressed in terms 
of the same universal scaling function g„ containing all the /U-dependence. 

Moreover, the results of this section clearly only depend on the parameters /i, u, T, so 
forgetting about the dependence of fx and u> on f3, we may in particular interpret the cases 
fx = 1/2 and li = 3/2 in terms of a one-dimensional harmonic oscillator. Indeed, when 
fx = 1/2 (resp. fx = 3/2), our quantum system reduces to the even (resp. odd) sector of an 
ordinary harmonic oscillator with frequency uona time circle of length T. More precisely, 
the eigenvalues E n = u(2n + fx) (|5.16|) and eigenvectors t(j n (|5.12|) reduce respectively to 



l Z 2 n 2 a/(2?7,)! 

/*=|: E n = u(2n+l + \ Vn(*0 = ~^== (-) 1 H 2n+1 {^)e~ ^ , 
z l 2 n y'(2n +1): 7r 

(5.52) 

in terms of the Hermite polynomials H n {x). Note that f is still restricted to be positive, 
but we may relax this condition by remarking that the Hermite polynomials of even degree 
are even, while the odd degree ones are odd. So we may extend the range of <p to the whole 
real line, upon simply redefining the eigenvectors above ip n — > ip n / a/2, and considering only 
even observables, say. The corresponding partition functions read respectively 

g 2 e 2 
Z T i = — - and Z T 3 = tt^, (5.53) 

and their sum is nothing but the partition function of the harmonic oscillator 

_ ujT_ 

e 2 

Zt, osc = Z T> i + Z T 3 = - - ^_ ujT . (5.54) 

More generally, the result (|5.49| ) for the correlation functions of the operator e w ^ z can 
be used to derive the corresponding quantity for the harmonic oscillator. Indeed, the odd 
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and even sectors are orthogonal, and the full correlator reads 

k 



m=l 

= sinh ' i 



1 



^T,osc 



iZ T igi{s) +Z T sgs(s)) 



= sinh 



2 j y i-s 2 

uT\ I 2s 



( S (l ~ Vl-S 2 )-^ + (1 - Vl-S 2 )^) 



(5.55) 



1 - s 



where s = 2/P(Z\, Z k \ti, t^) and as in (|5.50|) . From ( |5.55| ) we learn that the cor- 
responding correlation function in the harmonic oscillator case also only depends on the 
arguments of the observables through the quantity s. In Appendix C, we compute directly 
the correlator ( |5.55|) within the framework of the harmonic oscillator and obtain an alter- 
native expression for the combined argument s, which leads to the following determinantal 
expression for the polynomial P 



„ i i i 2/<*>2\ det(M-Z), 

P(Z X , Z k \t x , ...,t k ) = 2(1 + 2sinh 2 ( — ) ^ } ). 



(5.56) 



where M and Z are the following k x k matrices 



M 



( 7i -xi 
-xi 72 -x 2 

-x 2 73 



\-Xk 







fZ x 
Z 2 



-£fc-i 7fc 



/ 







V o 



z 3 



\ 





o z k J 



(5.57) 



with Xi = l/(2sinh(wti)), and 7^ = 2xiXi-ismh(u>(ti + i = 1,2, ...,k, with the 

convention xq = x k , to = t k - 



6. Discussion 

6.1. Lorentzian vs. Euclidean triangulations 

With the present work we have added yet another case to the list of problems for 
which the combination of methods from quantum field theory and statistical mechanics 
has proven very powerful (see for instance p[], ||14j| ). 



We have seen that in a certain region of the coupling constant space of our model, 
more precisely for f3 < 1, it is possible to define a continuum limit. The coupling constant 
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(3 survives this limit in unrenormalized form and we get a one-parameter family of contin- 
uum models which are neither identical to the continuum description of pure Lorentzian 
triangulations nor to that of Euclidean ones. The models do, however, have some common 
features with the pure Lorentzian triangulation model: they have the same scaling of the 
time variable and thus the same fractal dimension dn = 2 of the continuum space-time 
manifolds. One could have hoped, bearing in mind our original aim at gaining a better un- 
derstanding of the renormalization idea of || , that one would get out Euclidean quantum 
gravity for (3 = 1. However, we are led to the conclusion that no continuum limit exists 
for (3 > 1. The class of outgrowths that we considered here is too restricted to allow the 
model to flow from the Lorentzian case to the Euclidean one. What happens instead as (3 
increases from zero to one is that the correlation length £ defined by (cf. eqn. ( |5.17| )) 

Zt ~ exp ( — —T ] , as T — > oo, (6.1) 



i.e. 

decreases from to zero. This means that our time slices effectively decouple from each 
other, not allowing us to interpret our model as a model of surfaces. This decoupling effect 
is of course due to the fact that the outgrowths which dominate the geometries more and 
more as (3 increases do not imply any interaction between different time slices. 

6.2. Effective integrable structure and Renormalization Group formulation 

In ref. , a two-parameter family of theories has already been introduced to describe 
Lorentzian triangulations with intrinsic curvature energy (the parameters being the weight 
per triangle and a curvature weight). These models were found to be integrable, in the 
sense that their transfer matrices taken at different values of the curvature weight and 
of the weight per triangle but along the same spectral curve commute with one another. 
However, all these models turned out to share the same continuum limit up to a redefinition 
of the space and time scales which absorbed the (irrelevant) curvature parameter. One of 
those models is precisely the pure case (3 = 0. 

Here we have a very different situation: the parameter f3 survives and leads to distinct 
continuous theories. On the other hand, the model we started from does not seem to be 
part of an integrable family (except for (3 = 0). Still, the very simple form of the energies 
( |5.16|) suggests that for each value of (3 < 1 there exists a two-parameter family of transfer 
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matrices which is integrable and leads to the same continuum limit as our original model. 
The corresponding kernels are readily constructed by use of the quadratic relation ( |5.39| ) . 
Upon taking 

z = A, cp = v 7 ", ¥t = Va 7 , (6.3) 



we may introduce the kernel 

Gn, A (a,a ) = 2^ — 7^=^ , — y= A 2 



1 - A y 1 — A 



(6.4) 



Vaa' 



where \i — 1 = /3/(l + /3 2 ) as before and the ?/4 are given by ( 5.12 ) with the substitution 
uj — > O. By construction, it is clear that 



/"OO 

/ dot"Ga,\{ot., a")Gfi t x'(a",a') = Ga,\\>{oi,oi') 
Jo 

hence all the transfer kernels commute at fixed O. Choosing 



(6.5) 



1 + A \l + (3 

and rescaling the a's as 



fi = fi ( A ) = T^4fT^) > ( 6 - 6 ) 



«=( TT ^J «, a =(l^2 ) «. ( 6 - 7 ) 



we get a one-parameter family of kernels 



Gx(a,a') = G aWtX (a,c/) = J^e"*^)/^ (2^-^) . (6.1 

The kernels (|6.8|) lead to the same continuum limit a — > as our original generalized 
Lorentzian triangulation model with parameters f3 and g# upon taking 



A = e - 2 -« = e -wx(!±g). 



(6.9) 



with the same a and A as in ( |4.24| ). 
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We may now reinterpret the composition formula (|6.5|) as a self-similarity property 
for the kernel G\. Indeed, taking A = A', we get 



/"OO 

V dada' / da"G\(a, a")G\{a" , a) = Vdada'Gfi 
Jo 



(l + AF™ e VT ' ^V^l + A) 

1 + A 2 M \ 1 + A 2 J 



= V da,da'G\2 (a, a'), 

(6.10) 

where we have performed the change of variables a = n^fp ^ and similarly on a', and 
as before a = ( \+p2 ) 2 a. The operation performed in (|6.10|) decomposes into two steps: 
(i) composition of G\ with itself, by integration over the intermediate variable a" (ii) 
rescaling of the variables a, a' — > aa', to yield the kernel G\2. This is nothing but a 
real-space functional renormalization group (RG) decimation procedure. Indeed, imagine 
we wish to evaluate a product of say 2 n identical kernels G\. We first apply ( |6.10| ) to the 
2 n ~ 1 consecutive pairs of kernels, and end up after rescaling with a product of 2 n_1 kernels 
G\2, with the renormalized value A 2 of the parameter. Iterating the process will lead to 
the final single kernel G X 2 n (ot n , a' n ), where 



Qn= q A 2„- M2 «n-i, a = a, (6.11) 



and similarly for a' n . This is easily solved as 



1 + A 2T ' 1-A . 



and similarly for a' n . This yields an alternative derivation of the loop- loop correlator ( |5.35| ) 
over a time lapse t = 2 n . Setting A = e~ 2auj , t = 2 n = T/a as usual and taking the initial 
boundary values 

ao = a = \TVW) T' " = a=2 liT^Jv (03) 

we get the iterated values 

a n = 2a; cotanh(c;T) f - — ^ 2 j L , Q n = 2u cotanh(cuT) f - — ^ 2 j L T , (6.14) 
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so that the loop-loop correlator reads 



l-/3 2 ^ 2 



Zt(L 0i Lt) = 2 ( - | J u cotanh(t<;T) x 



(6.15) 



1 + /3 2 

G e -2^T ^2u cotanh(t<;T) ^ - — ^ 2 J Lo, 2w cotanh(u;T) ^ - — ^ 2 J Lt 

which is nothing but (|5.40|) . 

The self-similar kernels under RG transformations have been studied and classified 
in ref. [jT^J in the context of 1+1-dimensional critical wetting, and correspond precisely 
to the form (|6.8|) above. Each universality class is entirely characterized by the value of 
— 1 = P / (1 + P 2 ) . Our original model for decorated Lorentzian triangulations provides 
therefore an explicit discrete and simple realization of these universality classes. 

6.3. Calogero vs random walks 

For f3 < 1 where a continuum limit exists our generalized Lorentzian triangulation 
model is equivalent to a one-dimensional Calogero model. This model is known to appear 
in many different physical problems. It would be interesting to understand whether our 
generalized Lorentzian triangulation model thus has an interpretation in terms of a com- 
pletely different physical system. One particular realization of such a different equivalent 
physical system has been found in || where the pure Lorentzian case (3 = was shown 
to be equivalent to a one-dimensional Random Walk (RW) confined in a segment of size t 
upon identifying the fugacity per triangle g with the fugacity per step of walk. The RW 
model is known to be expressible as a massive scalar field theory in one dimension with 
action 

S(<p) = \J du{<p'(u) 2 +uj 2 ip(uf), (6.16) 

which is nothing but the ordinary quantum theory of a harmonic oscillator. The quantity 
ip(u) 2 is interpreted as the density of RW at point it, which precisely translates into the 
quantity Af(u) ( [3.271) under the abovementioned equivalence, while u is the continuous 
counterpart of the fugacity per step of walk. 

Surprisingly enough, we have stumbled here on the one- dimensional Calogero model 
at A = 1, rather than the harmonic oscillator corresponding to the action ( |6.16| ) above. 
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Fig. 5: Random Walk picture of pure Lorentzian triangulations. A 
Lorentzian triangulation (a) of time size t and its representation as a RW 

(b) between two reflecting walls at distance t + 1 (to go from (b) to (a), sim- 
ply squeeze the triangles horizontally). The reflected walk can be unraveled 

(c) , each portion between two contacts being reflected or not (here we made 
for instance two reflections at positions indicated by the dashed vertical lines. 
We finally compactify the time direction onto a circle of perimeter 2(t + 1), 
ending with a RW (d) attached at the antipodal points and t+1. 

Note that due to the extra potential l/(8<£> 2 ), a crucial difference is that the range of ip is 
over IR_|_ instead of IR for the harmonic case. 

It is an interesting exercise to re-derive the pure Lorentzian triangulation results in 
the RW harmonic oscillator language. In the following we will concentrate on the loop-loop 
correlator Z(i,j,t) defined in fl5.18p . In ref. 0, it was shown that 

Z itj (t) = g l Z{i -l,j- l,t) = J2( 2 9) Sp (hJ,S;t + 1) (6.17) 

^ ®> S>0 

where P(i, j, S;t+1) is the probability for a walker on the discrete integer segment [0, t + 
1] C IN making steps of ±1 with probability 1/2 at any integer position < h < t + 1 and 
being reflected with probability 1 whenever it reaches the boundaries h = or t + 1, to go 
from the position h = to the position h = t + 1 in S steps, with exactly % visits at h = 
and j visits at h = t + 1 (see Fig. 5). The slightly different transfer matrix g l+ ^~ l ( l_ ^^ 2 ) 
used in is responsible for the the shifts by —1 of i and j, and for the factor g per time 
step. 
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The reflection condition at h = and t + 1 may be unraveled by associating to each 
walk on [0, £ + 1] the set of walks on 2 obtained by iteratively reflecting or not each portion 
between two contacts at the boundaries (see Fig. 5) These new walks now go from the 
position h = to some position h = (2k + l)(t + 1) with fceZ. The numbers of original 
contacts i and j with respectively the boundaries h = and h = t + 1 translate into the 
total numbers of contacts with respectively the lines h = 2k(t + 1) and h = (2k + l)(t + 1), 
G Z. By compactifying the target space TL onto the circle Z/(2(i + 1)2), we end up 
with walks from /i = 0to/i = t + lona circle of perimeter 2(t + 1) as shown in Fig. 5. 
In this final picture, i and j represent respectively the numbers of visits of the antipodal 
positions h = and h = t + 1. 

In the continuum limit t = T ' j 'a, g = |(1 — a 2 ^), the latter RW's are governed by an 
action of the form (|6.16|) with u G [0, 2T] and the cyclic boundary condition <p(u + 2T) = 



<p(u). In this language the rescaled numbers of visits ai and aj become respectively ^(O) 2 
and <p(T) 2 . Moreover, it is simpler to compute the rescaled loop-loop correlator ( |5.22j ) 
with (3 = 0, by weighting each visit at h = and t + 1 respectively by x and y and setting 
x = 1 — aX, y = 1 — aY as before. We end up with the equivalence 

Z (X Y) - (M0MT)e- x ^ 2 e -rHTf } 

Z T (X,Y) — ——^——^ (6.18) 

where (...) stands for the functional average using the abovementioned action for the har- 
monic oscillator on a circle of perimeter 2T. In eqn.( |6.18| ), the factors <p(0) and <p(T) ensure 
that the walk starts at and ends at T (the factor 2 relates to our choice of normalization 
of the action). The insertions of exponential terms account precisely for the weights x,y 
per visit at 0, t+1, and the denominator ensures the suppression of the contribution of 
disconnected walks. It is now a simple exercise in ID quantum mechanics to evaluate the 
rhs of ( 6.18 ) namely, using the harmonic oscillator Hamiltonian H osc = — \-^z + \oJ 2 ip 2 : 

Ho ^od WT d VT e-^ x ^TY) (V|e- T/ H^o>) 2 
Z T (X, Y) = — ^ ^J-, (6.19) 



to be compared with the Calogero formulation (|5.23|) . Using the diagonalization formula 
(C.4) of appendix C below, with z = e~ ujT ', we are left with a ratio of two two-dimensional 
Gaussian integrals 

Z T (X,Y) = Ju - ^ 2 _ = 2(M" 1 ) 12 (6.20) 
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with the matrix 



i X + w cotanhfTcu) — . Jf /rr - s , 
M = 2 w v smh(Tl T ) ^ x , 6.21 

1 -ifflfcy F + ^ cotanh(Tu;) ' 



leading finally to the previous result ( |5.33| ) , with oo = v A. 

In view of the above derivation, it is remarkable that the Calogero formulation of the 
problem automatically takes care of all the subtleties (boundary conditions, suppression 
of disconnected loops, etc ...) we encountered in the RW approach. In this regard, the 
Calogero formulation of pure Lorentzian triangulations is better adapted. 

6.4 ■ Calogero vs 2D harmonic oscillator 

The previous section is a manifestation of the well-known connection between the 
Calogero Hamiltonian at fj, = 1 and the harmonic oscillator. As we will see now, the 
Schwinger representation we have developed in this paper provides another direct and 
exact link between the pure model and a pair of (untwisted/twisted) harmonic oscillators. 
This connection holds even before taking the continuum limit. 

Let us start again from the /3 = 0, h = 1 counterpart of ( |3.7| ) in which we also have 
incorporated for convenience an entropic factor g per time slice, namely 

g*ZW(t) = 11(1^1°° d ^~ aa ) f[ge 9aAPr+ ^\ Po = Pt- (6.22) 

We are about to show that the coordinates a > and p e £ may be rephrased into 
polar coordinates of the real plane, pertaining themselves to two one-dimensional harmonic 
oscillators. Performing the change of variable 



/V = 1/— e iXr+1 , r = 0,l,...,t-l, 

«r (6.23) 

a r = R^, r = 0, 1, t, 

with cto = ctt and where we have deformed the contours so as to have all \r £ [0, 27r) , we 
get 

g'Z^it) = i[( ^ 2R s dR s e- R ^\ J[ ge 2 ^ R ^ cos ^\ (6.24) 

s =l \Jo 77 JO J r=1 
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Finally, let us perform the following changes of variables on the x's, according to the parity 
oft: 

r + r _i, r = 2,3,...,t 

e 1 + e t , r = l 

(6.25) 



t odd : Xr 



t even : Xr 



r + r _i, r = 2,3,...,i 
01 - t , r = 1 



with 6 r G [0, 27r), r = 1, 2, The extra sign in the case t even is ad-hoc to make the 
change of variables invertible. We can now interpret 

(x r = R r cos 9 r , y r = R r sin r ), (6.26) 

as two real coordinates of the plane IR 2 , so that (|6.24|) becomes 



8=1 ^ ' r=2 



(6.27) 

This last formula suggests to introduce the following transfer kernel: 

G g (x,y;x',y') = 9 - e -U- 2 +y 2 +-' 2 +y' 2 ) e ^'-yy') ? (6 . 28) 

7T 

in terms of which the partition function reads: 

t „ t 
g t Z ( ®{t) = Yl( / dx a dy 8 )G g (x t ,(-l) t - 1 yt-,x 1 ,y 1 )Y[Gg(x r -i,y r -i]X r ,y r ). (6.29) 

s=l r=2 

The transfer kernel ( |6.28|) factorizes into G g (x,y;x',y') = G g ~(x,x')G g ~(y,y'), where the 
two kernels G^ respectively correspond to the harmonic oscillator and its twisted version, 
namely 



G+(X,X') = ^ji e -*(- 2 +-' 2 -4 9 Xx') = ^^080^^080^/)^+^ 

G-{y lV ') = e -4(y 2 +t/' 2 +W) = J2(-i) n ^n c (y)r n sc (y')q nH 



(6.30) 



7T 

n>0 



where we have set 



(6.31) 



and the n-th eigenfunction of the harmonic oscillator reads 



C SC M = X l -j^Hni^e-W , (6-32) 
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where H n (x) are the Hermite polynomials. The factorization of the kernel implies that of 
the partition function 

g'Z^it) = Zl sc (t) Zl sc (t), (6.33) 
where, writing formally products of kernels like ordinary products, we have 

Z^(t) = Tr((G+Y) = Z™(t) = TrdG-yS 1 ^) = (6.34) 

where S has the kernel S(y,y') = S(y + y'). Substituting these into ( ]6.33| ) yields ( |3.19| ). 
When we take the continuum limit a — > ( |3.23| ), both Z+ c and Z° sc become functional 
integrals of a scalar free field with action ( |6.16| ), the only difference being the boundary 
condition, namely <p(u + T) = <p(u) for the oscillator and ip(u + T) = —<p(u) for its twisted 
version. Indeed, in the latter case, we may perform the change of variables y^ r — > —y2r 
for the even indices, while the y^r+i remain unchanged in order to recover the correct 
sign of the untwisted case, but for both parities of t we end up with a wrong sign for the 
term yiyt, translating in the continuum limit into the above twisted boundary condition. 
The occurrence of both (twisted and untwisted) sectors arises from the fact that our 
boundary condition simply imposes that the rescaled number of triangles M at u = and 
u = T are identical, namely <p{u + T) 2 = <£>(«) 2 , which decomposes into periodic and anti- 
periodic boundary conditions for the free field itself. This factorization property extends 
to correlations involving the observables £ s (z) ( |3.12|) , namely 



m=l m=l 

k 



(6.35) 



tt e + m \ / TT e m \ 



This explains in particular the factorized form of the one-point average ( |3.22| ). Note that 
the above steps cannot be retraced for /3 > 0, where no such free field formulation seems 
to exist. 

It is instructive to re-derive the loop-loop correlator along the same lines as above. 
We start from 

g^Z^\x,y,t) = J] I ft / da s e-°X{ g e 9a ^ + — \ ± = x , p t+1 = y. 

(6.36) 



37 



Let us change variables as before to a r = R 2 , r = 1,2, t + 1, p r = e lXr+1 R r+ i/ R r , r = 
1,2, ...,t, and Xr = &r + &r-i, t = 2, 3, ...,t + 1, where we choose 6>i arbitrarily. We end up 
with t + 1 two-dimensional integrals over (x r , y r ) = (R r cos r , R r sin 6 r ), r = 1,2,...,£+1, 
by adding for convenience an extra integral ^ J Q n d9i, since the result does not depend 
on the choice 9\. This yields 



g^zM(x,y,t) 
t+i 



s=l 



t+i 

n 

r=2 



2g(x r -ix r -y r ^ 1 y r ) 



Note that if we switch to ?/2r —* V2r an d 2/2 
the y's are identical, leading to 



r-l 



(6.37) 

-j/2r-i, the integrals over the x's and 



^+ 1 zW(x, 2/ ,t) = det- 





-1 

-1 




\ 



-1 i- 



(6.38) 



VJ 



q\l-q 2 ) 



(1 - qx)(l - qy) - q 2t (q - x)(q - y) ' 

where - = q + |, and the matrix is of size (t + 1) x (t + 1). This is in perfect agreement 
with the result of ref. and with the continuum limit (|5.33| ). 

To conclude this section, we have unearthed an explicit two-component free field 
structure in the pure Lorentzian triangulation model, as opposed to the one-component 
free field description of Sect. 6.3 above. This increase in dimension is apparently the price 
to pay for reducing the interval [0, 2T] of the latter description to the original one [0, T]. 
Indeed the information needed to reconstruct the observables at time u pertains to the 
pair of slices u and 2T — u in the one-dimensional free field model. 

This two-dimensional structure a posteriori explains the emergence of the A = 1 
Calogero model in its continuum description: indeed, the corresponding two-dimensional 
Schrodinger equation 



(6.39) 
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Upon redefining ip(R\t) = R 2 0(i?|t), we end up with the Schrddinger equation with 
Calogero potential 




In this section, we have only considered radial observables thus ignoring completely the 
angular sector of the theory: it would be interesting to find the meaning of the angles in 
the triangulation language. 

Finally, it would be interesting to understand if it is possible to construct a surface (or 
volume) model for which the continuum theory would belong to the class of higher dimen- 
sional (multi-particle) Calogero models which are known to have a much richer structure 
than the simple one, but this is another story. 
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Appendix A. An alternative analysis of the transfer matrix 

In this section we show how to recover the Calogero Hamiltonian with A = 1 in the 
pure case using the strategy of ref. 0. We also explain why the same strategy becomes 
much more involved for f3 7^ and cannot apply to (3 > 1. 

The starting point is the composition law || 




(A.l) 



where Q^(x,y,t) = (©W) (x,y) with Q^\x,y) given by (|275|). Here we can choose the 
contour C to encircle the singularities of {x, which consist of a simple pole, and 

not those of 8^ (u, y, t). Inserting the scaling relations 
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and assuming the time evolution to be described in the scaling limit by a Hamiltonian H, 
one gets for the evolution on an elementary time lapse T = a: 



[l-aH + 0(a 2 ))^ = \i r^e« fx, -) ^) 

2 J 27TIL) \ IV J 



(A.3) 

where c is chosen so that the integration contour lies to the right of the pole at Z = X , 
and we have denoted by the same letter i/j(u> = 1 — aZ) = ip(Z). This gives 

H=(X 2 -A)^+X. (A.4) 

Writing tp(X) = J Q dLe~ x L tp (L) , L has the interpretation of the continuum loop length 
(cf. eqn. (|2.3|) ) and the inverse- Laplace transformed Hamiltonian H acting on functions 
of L takes the form 

6 = - L & -m + Ai - (A ' 5) 

with a flat measure on L for the wave functions tp(L). Then performing the simultaneous 
change of variables and functions 

l =\p 2 ^ #^ = vWy), (A.6) 
guaranteeing a flat measure on (p for (j)((p), we finally get 

which coincides exactly with the result ( |5.10| ) for (3 = 0. 

Now, let us try to apply the same strategy for (3^0. Evidently, we need first to 
evaluate the transfer matrix. For the present calculation it proves convenient to consider 
its symmetrized version, i.e. 

e(x,,)= / ^/^(eW)^f x ,J_)e( 2 )( Wl , W2 )(eW)^fJ_,,V (A . 8 ) 

Since the transfer matrix of the pure model has been explicitly diagonalized 0, it is 
straightforward to write down an expression for (0( 1 ^) 1 / 2 (x, y) 

(0 (1) ) 1/2 (x, y) = 1 = 1 Cl r iy \ , (A.9) 

1-C 2 (x)y 1-C 2 {y)x 
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where 

(l + 2g) 1 / 2 g(l + x) 

Ci x = r, C 2 m = 7- r. (A. 10) 

l + g(l-xY l + g{l-x) 

Here we have chosen to write (0 < ' 1 - ) ) 1/ ' 2 (x, y) in a form which explicitly exposes its singular- 
ity structure. Notice that (0 ( - 1 - ) ) 1 / 2 (x, y) is symmetric in x and y as it should be. Inserting 
the expression ( |A.9| ) for (G^ 1 )) 1 / 2 and picking the pole u)\ = C2(x) in the first term and 
u>2 = C2(y) in the last one, we end up with 

e(x,y) = C 1 (x)C 1 (y)e^(C 2 (x),C 2 (y)) 

CiWdjy) 1 

y/l - 46> 2 (C 2 (x))Vl - 49 2 (C 2 (y)) 2 1 _ A^c 2 { x )c 2 { y ) 

(A.11) 

From this expression for the transfer matrix we can read off various facts about the critical 
properties of our model. First, we see that there are several ways in which the transfer 
matrix can become singular. One possibility is that the argument of the square root 
vanishes. Assuming that as usual the critical valuesil of x and y are x c = y c = 1 this 
situation occurs when 

400=1. (A.12) 



Another possibility is that the denominator of the second factor in ( |A.11|) vanishes. This 
corresponds to the situation 

2gd(p+%)=l, and < 1. (A.13) 



The critical point (|A.13| ) is the one which reduces to the pure Lorentzian triangulation 
critical point when j3 — > (more precisely when h = 1 and 9 — » 0). The singularity given 
by ( |A.12| ) corresponds to the critical point of an isolated arch system. As long as (3 < 1 we 
always reach the critical point of (|A.13j ) before that of ( |A.12| ). For /3 = 1 the two points 



coincide and finally for /3 > 1 only the singularity (|A.12| ) persists. This analysis matches 
exactly our results in Sect. 4, in particular the relations (|4.14|) , (^4.17| ) and ( |4.23| ). 



5 The reader should not be surprised that the critical values of x and y differ from those of 
( 5.20| ). Indeed, the critical values w c = z c = of ( |5.20| ) correspond to the vanishing of the 



denominator of 6 (1) (>, z) = § ^ ^(® W ) 1/2 (w, y) r J ^(9 (1) ) 1/2 (x, z), while the denominator 
of (Q^) 1//2 (x, z) vanishes for x = x c = 1 when z = z c , and similarly for y = y c = 1 when w = w c 
(by use of flPo|) ). 
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It is obvious that the scaling behavior of 0(x, y) in the vicinity of a singularity cor- 
responding to the vanishing of the square root is not compatible with the existence of a 
Hamiltonian (cf. relation ( |A.3| )). Thus we can only hope to be able to introduce a mean- 
ingful continuum time variable as long as (3 < 1. This conclusion is of course in agreement 
with the analysis of Sect. 4. It is, however, not straightforward to write down the Hamilto- 
nian for (3 < 1 using the relations ( |A.1| ) and (|A.3| ) , the reason being that for the generalized 
model we will get not only a pole contribution to the contour integral in ( |A.3| ) but also a 
cut contribution. Introducing the scaling relations 

£=~fl-~a 2 A), x = l-aX, 6 = (b + , (A.14) 

the contribution from the pole in ( |A.11| ) reproduces the the same structure as that of ( |A.3| ), 
up to a redefinition of the renormalized fugacity per triangle, A -> A = A (jzf) ■ The 
contribution from the cut is more involved and cannot be obtained so simply. This is not 



surprising: the Hamiltonian (|5.10| ) (using (|5.8Q) should read in the language of modified 
lengths L=y 2 = (j^) L: 

H = -L—z =r +u 2 L + vr _ ' . (A.15) 

8L 2 dL AL 

When Laplace-transformed in terms of X, it becomes 

H V»(X) = ((X 2 - u 2 ) A + lV(i) + J™ ^( u )du, (A.16) 
showing that the cut contribution must give rise to a non-local integral term. 



Appendix B. The /c-point function of the scaled numbers of triangles for f3 < 1 

In terms of the time intervals t m = it m +i — u m for m = 1,2, ...k — 1, and tk = 
T — Yli< m <k-i tne ^-P°i n t correlator of e z muN(u m ) reac j s 

k k 

2 



(B.l) 

T,u Jo • i 
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where we have defined <f>k+i = fi- We now use k times the quadratic relation for the 
modified Laguerre polynomials ( p.38|) and the change of variables Vi = uxpf, i = 1, 2, k 
to rewrite 

( -Q e ^^„ l)) = 2sinh(u;T)e .T( M -i) / YlxidVi^e'^I^Xiy^v-^), (B.2) 

m=l ^° i=l i=l 

where we have set ^ 

2smh(wti) ( B3 ) 
7i = 2x^-1 sinh(w(ti + tj-i)), 

with the convention that to = tk- Let us evaluate by induction the integral 

k k 



J k {xi, x k \8 



r 

.,4) =2sinh(^T)e^- 1 ) / JJ Xl dv t \[ e~ 5 ^ I M _ 1 (2^V^m)- 

■ /o i=i i=i 

(B.4) 

We first perform the integration over v k by applying the following integral relation 

r ^ d«e- r %_ 1 (2 S ^Vi(2^) = e' (s2+s ' 2) /,-i J , (B.5) 

with r = <5fc, s = Xk-iy/Vk-i and s' = Xky/vi, with the result 

J fc (xi, scfc|5i, 8 k ) = J fc _i(xi, x 2 , — = \8 ± - 8 2 , 4-2, 4-i t~^)- ( B - 6 ) 

0fc Ofc Ok 

This is valid for all k > 3. When fc = 2, we simply get a single integral 
J 2 (xi,x 2 |(5i,(52) =sinh( W T)e^- 1 )^^ / dve~^~^ )v I^ 1 (2vx 1 x 2 /8 2 ), (B.7) 

^2 Jo 



readily evaluated as a particular case of the Laplace transform of the modified Bessel 

i± 

2x±X 2 



2 , 2 

function (|5.45|) with r = 81 — X2 and s = 2x1X2/82- This gives 



Ja(a;i, x 2 \8i,8 2 ) = 



8i8 2 — x\ — x\ ' 



9,(3) = sinh^rje^-V-^ 1 — 
Similarly, for /c = 1 we get 



(B.8) 



2ri 

= /„(*) = J l{xM = 9 ,(^) = g ,( cosh( , r) _^ sinh( . T) ), (B.9) 
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In general, applying the recursion relation ( [B.6Q , we find 



( fi ^*. )=9 , { J ), (B.!0) 

m=l 

where the polynomial P is defined as follows: we first introduce the sequence 

x^ 

fJ>k- P = Sk-p p = l,2,—,k-l with fJ,k = 8k, (B- 11 ) 

Mfc-p+i 

and the polynomial P reads 

k-2 

T>(<7 rz rz U + + \ V-XV-l—Vk «2«3---Mfc-p X k -pXk-p+l---X k , , 
ry&X-, £>2, ■■■,^k\tl,t2, ...,tk) = > , X . 

xix 2 ...a;fc xix 2 ...Xfc-p-i M fc -pMfc-p+i-"Mfc 

Once expressed in terms of and t^, this polynomial takes the very simple form ( |5.51|) , 
with for instance 

k = 0: P = 2cosh(^T), 

k = 1 : PiZilh = T) = 2(cosh(wP) - Z^inhfuT)), 

(B.13) 

k = 2 : P(Z U Z 2 \ti, t 2 = T-ti) = 2(cosh(wT) - (Z x + Z 2 )sinh(wP) 
+ 2i?ii? 2 sinh(i i ;ti)sinh(i ( ;£ 2 )) , 
leading respectively to the correlators 

k = 0: (1)=^(1) = / M (0) = 1, 

* = 1= (^)=^(p(|^)=/M(2), 



/„ Z 1 +Z 2 -2Z 1 Z : 



P(Zi, Z 2 |w 2 - «1, P + Ml - li 2 ) 

sinh(ci;(w 2 — Wi))sinh(c<;(P + ui — w 2 )) 



sinh(wP) 



(B.14) 



with as in ( |5.46| ). More generally, the /c-point correlator is expressed as 



2, P fc 



^ (« = ^) = /, (cotanh(u;P) - — _) , (B.15) 

in terms of P fc = P(i?i, ...,Z k \tx, ...,t k ). 

As a final check on the result ( |5.51|) , we see that 

P(Zi, Zk-i, Z k = 0|£i, ...,£fc) = P(^i, Zk-i\ti, tk-i) 
P(Zi, ...,Zk\ti, ...,tk-i, tk = 0) = P(Zi, Z k -2, Zk-i + 2k\ti, tfc-i), 
as expected from the definition of the correlator. 
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;b.i6) 



Appendix C. Comparison with the harmonic oscillator 

We wish to compute the /c-point correlation function of the operator e Uip ^ z in the 
one-dimensional quantum system with Hamiltonian 

on a time circle of length T. The Hamiltonian has eigenvalues E n = u(n + i) and eigen- 
functions 



where H n (x) are the Hermite polynomials, n = 0, 1,2, ... and 99 G (—00, +00). Next we 
will use the following quadratic relation obeyed by the Hermite polynomials 

~ (z/2) n 1 2 X yz-z*(x2+y2) 

Y,H n {x)H n {y) K -^L- = -T^e I=* , (C.3) 

n=0 ' ^ 

immediately translated into 



£ vcwvr^i)^ = /_ ^-e*^^-T^*^+^). (a4) 



v ^(l-Z 2 ) 

n>0 v V ' 

The correlation function reads 

k 

( J^J e <+><p 2 (u m )Z m ^ osc 
m=l 

1 f fc 



r,osc J]R ^„ . , 



(C.5) 



C c (^)C c (^+i)^ 



where the time intervals £j = — itj, i = 1,2, fc — 1, ife = T — Xli<i<fc-i Applying fc 
times the relation ( |C4| ) with respectively <po, (pi — > (pi, (pi+i and z — ► e - ^, and performing 
the change of variables = <£>j, we get 

(f[ e VK)^) o8c= Vsi-Sfc /" (ff-^L) e -IE- =1 ^^ 2 - 2 ^^^, (C.6) 

mil -/R tl^ 

with Xi and 5j as in ( |B.3| ). This finally yields 

( II e^ 2 ^ z -) osc = 2sinh( W T/2) v ^TZ^det(M - Z)"* = ^ z y ( C - 7 ) 

as a straightforward result of the Gaussian integration, with M and Z as in ( |5.57|) . The 
comparison with the alternative expression ( |5.55|) immediately yields the determinantal 
expression for the polynomial P ( |5.56| ). 



45 



References 



[1] J. Ambj0rn, B. Durhuus and T. Johnsson, Quantum Geometry, Cambridge University 
Press, 1997. 

[2] F. David, Nucl. Phys. B257 (1985) 45; V.A. Kazakov, Phys. Lett. 150B (1986) 140. 

[3] J. Ambj0rn and R. Loll, Nucl. Phys. B536 (1998) 407, |hep-tli/9805108| . 

[4] R. Nakayama, Phys. Lett. B 325 (1994) 347. 

[5] B. Duplantier and I. Rostov, Phys. Rev. Lett. 61 (1988) 559. 

[6] J. Ambj0rn, J. Correia, C. Kristjansen and R. Loll, Phys. Lett. B475 (2000) 24, fiiep- 
| th/9912267| . 

[7] B. Duplantier and I. Rostov, Nucl. Phys. B340 (1990) 491; B. Eynard, E. Guitter and 
C. Rristjansen, Nucl. Phys. B528 (1998) 523, |cond-mat/9801281| ; S. Higuchi, Mod. 
Phys. Lett. A 13 (1998) 727, |cond-mat/9806349| . 



[8] J. Ambj0rn, R.N. Anagnostopoulos and R. Loll, Phys. Rev. D60 (1999) 104035, |hej> 

| th/ 9904012] ; Phys. Rev. D61 (2000) 044010, |hep-lat/9909T29[ 

[9] P. Di Francesco, E. Guitter and C. Rristjansen, Nucl. Phys. B567 (2000) 515, |Eep^ 

| th/9907084 . 



[10] V.G. Rnizhnik, A.M. Polyakov and A.B. Zamolodchikov, Mod. Phys. Lett. A3 (1988) 
819. 

[11] L.D. Landau and E.M. Lifshitz, Quantum Mechanics, Pergamon Press, 1958. 
[12] H. Bateman, Higher Transcendental Functions, Vol. I, McGraw-Hill (1953). 
[13] H. Bateman, Higher Transcendental Functions, Vol. II, McGraw-Hill (1953). 
[14] C. Itzykson and J.-M. Drouffe, Statistical Field Theory, Cambridge University Press 
1989. 

[15] H. Spohn, Fixed points of a functional renormalization group for critical wetting, Eu- 
rophys. Lett. 14 (1991) 689. 



46 



